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ABSTRACT 



Aims. With new THz maps that cover an area of ~3. 3x2.1 pc^ we probe the spatial distribution and association of the ionized, neutral 
and molecular gas components in the Ml 7 SW nebula. 

Methods. We used the dual band receiver GREAT on board the SOFIA airborne telescope to obtain a 5' .7x3'. 7 map of the '^CO 7=13- 
12 transition and the [C ll] 158 /im fine-structure line in M17 SW and compare the spectroscopically resolved maps with corresponding 
ground-based data for low- and mid- 7 CO and [C l] emission. 

Results. For the first time SOFIA/GREAT allow us to compare velocity-resolved [C ll] emission maps with molecular tracers. We see 
a large part of the [C ll] emission, both spatially and in velocity, that is completely non-associated with the other tracers of photon- 
dominated regions (PDR). Only particular narrow channel maps of the velocity-resolved [C ll] spectra show a correlation between the 
difi"erent gas components, which is not seen at all in the integrated intensity maps. These show difi"erent morphology in all lines but 
give hardly any information on the origin of the emission. The [C ll] 158 yum emission extends for more than 2 pc into the M17 SW 
molecular cloud and its line profile covers a broader velocity range than the '^CO i= 13-12 and [C l] emissions, which we interpret as 
several clumps and layers of ionized carbon gas within the telescope beam. The high-i '^CO emission emerges from a dense region 
between the ionized and neutral carbon emissions, indicating the presence of high-density clumps that allow the fast formation of hot 
CO in the irradiated complex structure of M17 SW. The [C ll] observations in the southern PDR cannot be explained with stratified 
nor clumpy PDR models. 
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1. Introduction 

Ml 7 SW is a gian t molecular cloud at a distance of -1.98 
kpc (IXu et al.ll201lh . illuminated by a highly obscured (A,, > 
10 mag) cluster of several (among ^ 100 stars) O/B stars 
(iBeetz et alJ 119761: lHansonetaT]ll997h. and it harb ors several 
candidate young stellar objects (iPovich et al.ll2009h . The large 

■ amount of observational data available in the literature and its 
nearly edge-on geometry make Ml 7 SW one of the best-studied 
prototypes of a clumpy photon-dominated region PDR inter- 

] face in the Galaxy. Studies of molecular and atomic emission 
indicate that the structure of the gas is highly clumped and 
not homogene ous (Stutzki e7al.lll988t IStiitzki & Gustenlll990t 
iMeixner et al.il 19921: iPerez-B eaupuits et al.ll2010i and references 
therein), and the structure of its neutral and molecular gas seems 
to be dominated by magnetic rather than b y thermal gas pressure , 
in contrast to many other PDR regions dPellegrini et al.ll2007l) . 
Temperature s of ~275 K were found t oward the VLA 21cm con- 
tinuum arc ( Brogan & Trola nJ'2001). and are associated with 
NHj and highly e xcited '^CO emission (Gtisten & Fiebig 1988; 

lHarrisetal.ill987l) . ^_^ 

Recent AKARI observations dOkada et alJ201(]|) with spatial 
resolutions between 39" and 57" showed that the [C ii] 158 /im 



line emission is widespread in the M17 complex and peaks at the 
northern (Ml 7 N) and southern (Ml 7 SW) bars, in agre ement 
with previous observations (Matsuhar a et al.|[l9 89: Stutzk iet al.l 
119881) . However, this and previous observations of the hot and 
ionized gas in M17 SW are limited in spatial resolution and 



exteiit (e.g. [Harris et al.l 1 19871: Istutzki et al.ll 19881: iGenzel et al.1 
1 19881: IMeixner et al.lll992l:lHowe et alj|200 0'). Therefore, in this 
work we present new observations (of an area 5'. 7x3'. 7) of hot 
molecular ('^CO 7=13-12) and ionized atomic ([C ii] 158 //m) 
gas, with spatial resolutions of ~19".8 and ~15".6, respectively, 
which advances existing work in M17 SW. 



2. Observations 

The observations were performed with the German Receiver for 
Astro nomy at Teraheitz Frequencies (GREAlQ, Hevminck et aLJ 
1201 2|) on board the Stratospheric Observatory For Infrared 
Astronomy (SOFIA). We used the dual-color spectrometer 
during its first Short Science flight on 05 April 2011 to 
simultaneously map the fine-structure transition of [C ii] at 
1900.536900 GHz (157.7yum) and the I^CO 7=13-12 transition 
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' GREAT is a development by the MPI fiir Radioastronomie and 
the KOSMA/ Universitat zu Koln, in cooperation with the MPI fiir 
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Fig. 1 : Left panels - Intensity maps of [C ll] 158 //m (top) and '"CO 7= 13-12 (bottom) 
in M17 SW, integrated in the velocity ranges 0-40 kms"' and 14-28 kms"', respec- 
tively. The contour levels are the 10%, 25%, 50% (thick line), 75% and 90% of the peak 
emission. Dashed lines are the strip lines shown in Fig.|2] The central position (0",0") 
is marked with a cro ss. The ultracompact H ll region M17-UC1 and four H2O masers 
JJohnson et alJll99a) are marked by the bla ck circle and plus symbo ls, respectively. 
Right panel - 21 cm continuum emission bv lBrogan & Trolandl ( 120011) with the over- 
laid contours of the velocity integrated (same as above) emission of '^CO 7=13-12 
(gre en), [C n] (red), and the [C l]^ P2 - ^^1 370 /im (blue, integrated in 9-30 km s"') 
from'Perez-Beaupu its et al.lll2010h . The contour levels (from thin to thick) are the 50%, 
75% a nd 90% of t he peak emission. The stars indicate the O and B ionizing stars 
dBeetz et al.lll976l : IHanson et alJll997r) . Dashed frames depict the beam center for the 
edges of the 6 OTF strips. Contour maps are smoothed to 20" resolution. 
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Fig. 2: Left - Strip lines of the molecular ('^CO) and atomic ([Cl], 
[C II]) gas at AS = 30" (PA. 90°) across the ionization front of 
M17 SW. Right - Strip lines at PA. 63° (AS = Att/2). The lower- 
7 lines of '-CO, and [C l] are from iPerez-Beaupuits et al.l ( 120101 
their Fig. 3). All these profiles are in units of K km s~' . The VLA 
21 cm continu um and H I optical depth ( integrated between and 
30 kms"') by iBrogan & Trolandl (1200 ih are in units of 500xJy 
beam"' and 8xt(HI) kms"', respectively. The offset, Aa = 0" in 
R.A., is the same as in Fig.[T] 



at 1496.922909 GHz (200.3 fim) toward M17 SW. The obser- 
vations were performed in on-the-fly (OTF) total power mode. 
The area mapped consists of six strips, each covering 224" x 32" 
(AaxAJ with a sampling of 8" (half the beamwidth at 1.9 THz). 
Hence, each strip consists of four OTF lines containing 28 points 
each. We integrated Is per dump and 5s for the off-source refer- 
ence. 

All our maps are centered on R.A(J2000)=18:20:27.6 and 
Dec(J2000)=-16:12:00.9, which corresponds to the SAO star 
161357. For better system stability we used a nearby refer- 
ence position at offset (345", -230"). A pointed observation of 
this reference position agai nst the k nown better reference (off- 
set (1040",-535") Matsuha ra"etan [T989) showed that the ref- 
erence is clean of '^CO emission, but contains weak (<20% of 
the peak emission) and relatively narrow (FWHM~15 kms"') 
[C 11] emission. All [C 11] spectra presented here were corrected 
for this missing flux. 

Pointing was established with the SOFIA optical guide cam- 
eras, and was accurate to better than 10". Because t he acousto- 
optica l and fast Fourier transform spectrometers (iKlein et al.l 
I2OI2I) operated in parallel give redundant information, we base 
the following analysis on the data from the latter, which pro- 
vided 1.5 GHz bandwidth with about 212 kHz (-0.03 km s ') of 
spectral resolution. Because during this commissioning flight the 
instrument showed some random gain variations, we discarded 



15-20% of the spectra by filtering the spectra for outliers in the 
total-power IF level and for outliers in the noise rms (obtained 
after subtracting a third-order polynomial), retaining only reli- 
able data. The calibration of these data to antenna temperature 
was perf ormed with the kal ibrate task from the kosma software 
package dGuanet al.ll2012h . We then reduced and imaged the 
data further with the CLASS90 packages, which is part of the 
GILDASQ software. Using the beam efficiencies (77^) 0.51 for 
[C II] and 0.54 for '^CO 7^13-12, and the forward efficiency 
(rjf) of 0.95 ( Hev minck et al J |20 1 2h . we converted all data to 
main beam brightness temperature scale, Tg - rjf x T^/rje- 

3. Results 

The velocity-resolved spectra of the new GREAT/SOFIA ob- 
servations reveal a much more complex structure than expected 
based on our understanding of Ml 7 SW as a highly clumped 
PDR, with a significant fraction of the [C 11] emission not at all 
associated with the molecular PDR material. 

Figures 1-4 display the GREAT data compared with supple- 
mentary data of typical PDR tracers from the literature. Figure[T] 
shows the velocity-integrated emission of [C 11] between and 
40 kms"' (peak -1070 Kkms"') and '^CO 7=13-12 [14 - 
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28 kms ', peak -180 Kkms'] (left panels). While the high- 
excitation CO emission follows the hot dense gas structure seen 
in previous submm line studies (but see the discussion below), 
the [C ii] distribution reveals surprises: its velocity-integrated 
distribution is quite shallow and far more extended than pre- 
dicted for the stratification in classical PDR models. This is 
demonstrated in Fig. [Ta] and the intensity cuts across the PDR 
in Fig. |2] showing strong [C ii] emission over the whole extent 
of the mapped area. 

The broader distribution is cons istent with th e lower angu- 
lar resolution data from the KAO dStutzki et al.| [T988), at the 
time interpreted to result from the deeper UV-penetration into 
a clumpy medium. However, the new GREAT/SOFIA data, re- 
solving the spectra in velocity, show a more complex scenario. 
We find that (in M17 SW) the [C ii] emission (£,,=91.21 K, and 
Hcrit ~ 5 X 10^^ cm"-') traces - in addition to the dense PDR gas 
west of the ionization front - also a widespread (toward the east), 
probably more diffuse atomic gas component that is invisible in 
CO. A sample of this is shown in Fig. [3] (left panel), display- 
ing in gray scale the intensity distribution of the lower velocity 
components of [C ii] and its (non)correlation with CO and H i. 
Comparing individual spectra (Fig. |4), we see that the span of 
observed velocities in [C ii] is much wider than that of any other 
PDR tracing species. In the spectra at offset (0,0), a position well 
in front of the PDR interface (Fig. [TbJ, for example, [C ii] still 
shows a bright and amazingly wide line with AV ~ 30 kms"', 
including velocities not detected in molecular (CO) or (denser) 
atomic (neutral carbon [C i]) gas column density tracers. Most 
interestingly, those velocities (e.g., 4 -10 km s~^) exhibit appre - 
ciable opacities in atomic hydrogen (Brogan & Troland 2001 - 
note that H i is observed in absorption against the 21 cm back- 
ground continuum, hence its detectability across the map de- 
pends on the continuum distribution). The broad and intricate 
structure of the [C ii] spectral line (Fig.|4]i can also be the signa- 
ture of several clumps or layers of ionized carbon gas at "redder" 
and "bluer" velocities than the other tracers. 

Fig.[Tb]provides evidence for a C^ — > C stratification across 
the PDR: the bulk of the [C i] emission (defined as >75% of 
the peak emission) avoids the bulk emission of [C ii], with 
their peaks being separated by ~70" in the N E-SW direction , 
i.e., 0.67 pc at the distance of -1.98 kpc (IXu et alJ 12017) . 
In a previous analysis the extended [C i] emission (west from 
the ionization front), and its particular peak emission in front 
of the bulk CO emission (if stratified, the [C i] peak should 
be observed before the CO peak, from the direction of the 
PDR interface), was argued to emerge from the interclump re- 
gions of a very clumpy medium, but could also result from 
a partial face-on illumination of the molecular clouds (e.g. . 



Stutzki et alJl l988l:IStutzki & Gusten|[T990llMeixner et al.lll992t 
Pere z-Beaup uits et al.l 120101) . In the context of the dumpiness 
of the PDR and the newly discovered association of much of 
the [C ii] emission with atomic gas (see below) at velocities not 
matched by any molecular material, the [C ii] peak cannot be 
interpreted as the classical stratification expected in a homoge- 
neous PDRQ but instead has to be caused by excitation gradi- 
ents with increasing optical depth into the atomic and molecular 
cloud. 

The excitation study of the warm PDR layers, including the 
CO excitation, will be the subject of a forthcoming publica- 



'' In a clumpy PDR, with many internal surfaces, the [C ll] peak 
should be even closer to the [C l] peak, becoming indistinguishable 
from it depending on dumpiness of the source and spatial resolution 
of the observations. 



tion. Here we briefly elaborate on the new GREAT data: The 
'"CO 7=13-12 line follows a similar spatial distribution as the 
lower-7 lines reported in Perez-Beaupuits et al. (2010), although 
the whole '^CO 7=13-12 emission is shifted toward the ioniza- 
tion front. Its peak emission is ~16" (0.15 pc) southeast of the 
'^CO 7=6-5 peak. The 7=13-12 line is already fainter than the 
lower-7 transitions, indicating subthermal excitation. 



4. Discussion 

We can distinguish two different spatial regimes in the PDR 
with embedded star cores of active star formation. The strip 
line along AS = 30" (PA. 90°, Fig. |2]\) traces the embedded 
star- formation, going through the ultracompact H ii region UCl 
(Aa = -30"), the embedded H2O maser positions farther west, 
and the weaker additional H 11 region seen as a separate clump at 
Aa = -1 10" in the 21 cm continuum map and matching the sec- 
ondary peak of [C 11], [C i], and the '^CO 7=2-1 lines (Figs.fTbj 
and|2}\). A dominant internal heating is likely suppressing the 
layering here. The strip with the position angle 63° (Fig. |2ji) 
instead traces the classical PDR. For a homogeneous PDR one 
expects a stratifi ed layering o f [C 11I [C i], high-7 CO and low- 
7 CO. However. IStutzki et al] (119881) showed the structure to be 
very clumpy so that no stratification should be observed. Our 
observations now show that high-7 CO and low-7 CO peak at 
similar depths in the cloud while the integrated [C 11] peak is 
shifted by x; 0.19 pc relative to the CO peaks. The [C 11] shift 
cannot even be explained by homogeneous PDR models, much 
less by a clumpy medium. 

In [C 11] we detect all velocity c omponents that were dis- 
cussed in the H i channel maps by iBrogan & Trolandl (1200 ll 
their Figs. 4 and 5), including a weak foreground cloud at 
7 kms , shocked gas in front of the cluster that is blown to- 
ward the observer at 11-17 kms ', and the main PDR velocity 
of ~21 kms '. In the higher density gas (west of the H 11 re- 
gion), [C 11] seems to be optically thick, showing self-absorption 
dips at the peak velocities of H i and the molecular lines. The 
'^CO 7=2-1 emission also traces the optically thin shocked gas 
down to 10 kms ' (Fig.©, the '-CO 7=13-12 line is hardly ex- 
cited there. A good match between [C 11] and '^CO 7=2-1 can 
be seen at intermediate velocities (Fig.|3] right panel) where the 
shifted [C i] emission with respect to [C 11] is also evident. 

We can compare the column density contributions visible in 
the different velocity components using the column density for 
the H I gas by Brogan & Troland (2001) and the C^ column den- 
sity (e.g., Schneider et al.i i2003t their eq.A.4) responsible for the 
observed [C 11] emission. For the gas temperature and the density 
we assumed 250 K and 10'' cm^^ at the off^set (0,0), and 150 K 
and 10^ cm~^ for the offsets deeper in the cloud according to 
iPerez-Beaupuits et all (l2010h . 

To compare our results with the H i data we estimated [C 11] 
and H i column densities at selected positions along the south- 
ern PDR (Fig. nil from the [C 11] emission and H i optical depths, 
integrating the spectra in the velocity ranges 11-17 kms"' and 
17 -24.5 kms~'. Selecti on of these velocity ranges is described 
bv ' Brogan et al.l (Il999l) . In the range 11-17 kms ' we found 
A^([C 11]) (with uncertainties of -20%) between 6.5x10'^ cm"^ 
at offset position (0,0) and 1.6x10"^ cm"^ at (-60, -30), while 
the corresponding H i column densities vary (with uncertain- 
ties 7-16%) between 5.5x10^' cm"^ and 4.7x10^' cm"^, re- 
spectively. In the velocity range 17-24.5 kms"' the [C 11] col- 
umn densities are between 1.2x10'^ cm"^ and 2.9x10'^ cm"^ 
at the same offset positions, and we found A^(H i) between 
3.7x10^' cm"^ and 4.4x10^' cm"^. The variation in HI column 
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: Selected 1 kms" wide channel 
maps (gray scale) of the [C ll] velocity- 
integrated (Kkms"') in the velocity ranges 
9.5-10.5 kms"' and 18.5-19.5 kms"'. The 
contours show the integrated emission in the 
corresponding velocity channels of [C l], 
'-CO 7=2-1, and t(HI) in 20% steps (from 
thin to thick lines) of the peak channel in- 
tegrated values (top right in the maps). All 
maps but r(H l) were smoothed to a resolu- 
tion of 20" . 



densities is less than 20% in both velocity ranges, indicating a 
relatively homogeneous H i interface seen also by [C ii] . The lat- 
ter, instead, presents a much stronger column variation (40%- 
60%) depending on the ambient conditions. 

Assuming the nor mal cosmic carbon abu ndance of 
X(C/H) = 1.2 X 10^* dWakelam & HerbsJ l2008h and com- 
plete ionization of the carbon, we obtained gas column densities 
A^(H) = (5.4 - 13) X 10^' cm"^ for the low- velocity compo- 
nent and A^(H) = (10 - 24) x 10^' cm"^ for the high- velocity 
component. The [C ii] emission of the low-velocity component 
in front of the PDR is therefore consistent with pure H i gas, 
while for all other components (above 17 km s~^) at least half of 
the [C ii] emission stems from molecular gas, i.e., classic PDR 
material. The material in front of the southern PDR must be 
ato mic to a high degree, indicatin g non-equilibrium chemistry 
(e.g. lStoerzer & Hollenbachll998l) . 

5. Final remarks 

The integrated intensity maps do not tell us about the origin of 
the emission lines. Only narrow channel maps of the velocity- 
resolved [C ii] spectra show a correlation with the diffuse gas 
components. From the comparison of different tracers in differ- 
ent velocity components we found that for the regions of embed- 
ded star formation the internal heating is always the dominant 
process, providing a very good match of all tracers without sig- 
nificant layering. The good spatial correlation of CO, [C ii], and 
other tracers at particular velocity c omponents/chann e ls con - 
firms the clumpy PDR picture from IStutzki & Giistenl (Il990h . 
However, a significant fraction of the [C ii] velocity components 
is not associated (spatially) with the dense (high-7 CO) or the 
diffuse (H i and CO 7=2-1) PDR tracers. The structure neither 
matches a stratified nor a clumpy PDR model. The [C ii] emitting 
gas is only partially explained by atomic gas. Optically thin, hot 
molecular gas must contiibute more than half of the emission. 
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